Self-consistent field theory for the 

o : 

^ : distal ordering of adsorbed polymer: 

^ '. comparison with the Scheutjens— Fleer model 

o ■ 
(N : 

; Karl Isak Skau, Edgar M. Blokhuis, and Jan van Male 

^ , Colloid and Interface Science, Leiden Institute of Chemistry, 

^ . 

^3 ' Gorlaeus Laboratories, P. O. Box 9502, 2300 RA Leiden, the Netherlands. 

S: 

^ ' February 2, 2008 

O ■ 

y—i ' Abstract. In a previous article [E.M. Blokhuis, K.I. Skau, and J.B. Avalos, 
> ■ 

CN ' J. Chem. Phys. 119, 3483 (2003)], a self-consistent field formalism was derived 
^ ! 

. for weakly adsorbing polymers, valid for any chain length. It was shown that the 

Q ! presence of a solid wall induces an ordering of the polymers on the scale of the 

^ ' radius of gyration far away from the surface (the distal region). These oscilla- 

g ' tions in the polymer concentration profile were first noted in work by Semenov 

. et ai, and later observed in numerical solutions of the Scheutjens and Fleer self- 
C ■ 

Q I consistent field model. In the present paper, we compare the weak adsorption 
O ■ 

• model in more detail with the numerical results from the Scheutjens and Fleer 

^ ■ model. Quantitative agreement is obtained for the polymer segment density pro- 

H ; 

5t ! file in a good (athermal) solvent. For theta solvent and poor solvent conditions, it 
is necessary to extend the weak adsorption model to take the non-local character 
of the polymer-solvent interaction into account. Again, quantitative agreement is 
obtained for the polymer segment density profile, in particular for the transition 
from an oscillatory decaying profile to a monotonically decaying profile when the 
bulk polymer density is below a certain threshold value. 
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1 Introduction 



The adsorption of polymer onto a solid surface has received considerable attention 
both from an experimental and a theoretical point of view. Interest in these 
systems is driven by practical applications, but also as a testing ground for various 
theoretical approaches [1-3]. 

Theoretically, surfaces with enhanced polymer adsorption were first studied 
by de Gennes [1,4] in the context of the Edwards self-consistent field theory [5]. 
In the de Gennes free energy functional, the so-called ground state dominance 
approximation [1,5,6] is made in which the polymer chain length is essentially set 
to infinity. Various extensions of the dc Gennes model have been formulated. In 
work by de Gennes [4] and Rossi and Pincus [7] the correct scaling behavior was 
incorporated into the free energy functional, whereas Semenov and coworkers have 
extended the de Gennes model to determine finite chain length corrections [8-12]. 
The latter calculations are of particular interest since the polymer chain length 
is an important parameter in experiments [13] and computer simulations [14-16]. 
Variation of the chain length thus provides a more stringent testing of theoretical 
models [17]. 

A versatile model for the theoretical description of polymer adsorption is the 
Scheutjens and Fleer [18, 19] self-consistent field model. The Scheutj ens-Fleer 
model (SF) is a generalization of the Flory-Huggins mean-field lattice model [13] 
extended to describe inhomogeneous polymer systems. In the SF model, the 
Edwards [5] self-consistent field equations for the polymer's Green function are 
solved numerically on a lattice. The capabihty and fiexibihty of the SF model has 
been demonstrated with its use to calculate properties of inhomogeneous polymer 
systems containing homopolymers, copolymers and polyelectrolytes and also sur- 
factant systems like micelles and membranes [2] . In recent years, the SF algorithm 
has been made available in the multipurpose computer program SFBOX [20]. 

Recently, a self-consistent field formahsm was derived for weakly adsorbed 
polymers, valid for any chain length [21]. One notable difi^crcnce from the clas- 
sical ground state dominance result of de Gennes [1] is the distal ordering of the 
polymers; the finite chain length gives rise to oscillations on the scale of the radius 
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of gyration away from the surface. These oscillations in the polymer concentration 
profile were first noted in work by Semenov et al. [8, 11], and later observed [22] 
in the numerical solutions of the Schcutjcns and Fleer self-consistent field model. 
The onset of these oscillations has also appeared in the context of analytical work 
and computer simulations [15,21,23]. 

In the present paper, we compare the free energy functional formalism for 
weak adsorption [21] with the numerical results from the Scheutjens and Fleer 
model. In particular we address the observation made by van der Gucht et al. [22] 
that in a theta solvent, the oscillatory decaying function becomes a monotonically 
decaying function when the bulk polymer density is below a certain threshold 
value, similar to the Fisher- Widom transition [24, 25] in the damped oscillations 
found in a simple liquid near a surface. Our comparison will therefore consider 
some general properties of the distal ordering in both good and theta solvent 
conditions. 

The outhne of this article is as follows: In the next two sections we discuss 
the self-consistent field theory for weak adsorption, in the context of which our 
calculations are made, and the Scheutjens and Fleer self-consistent field model. In 
Section 4, we make the comparison between our analytical results and numerical 
calculations carried out with the help of the computer program SFBOX. We end 
with a discussion of results. 



2 Self-consistent field theory 

The Green function G(r, r', N) describes the statistical weight of a single polymer 
chain of length N with one end at r and the other end at r'. The Green function 
is determined by the Edwards equation [5] : 

^ ^(r, r n) = ^V^ G(r, r', n) G{r, r', n) , (2.1) 
on 5 /cb-£ 

where a is the polymer segment length and C/(r) is an as yet unspecified external 
potential. As an initial condition to the Edwards equation, we have that 

lim G'(r,r',n) =(5(r-r'). (2.2) 
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In terms of the Green function, one can construct the average segment density: 



!^dn dr' dr"G(r',r,n)G(r,r'\N -n) , ^ 

(/. r = ATp ^-^ i ± c \ , 2.3 

where A^p is the total number of polymer chains. This prefactor determines the 
scale of the Green function; it is chosen such that the density is equal to the given, 
uniform bulk density, 0(r) = a^N^N/V = (pb, for a homogeneous system. Here, 
(f){r) is the polymer segment number density made dimensionless by the factor 
and can thus be interpreted as the polymer segment volume fraction. 

The statistical weight G{r, n) of a polymer chain with one end at position r is 
obtained by integrating the Green function over one of the chain's ends: 

G,...)./d.'G,.y,.). ,2.4, 

For the case of polymer adsorption onto a planar, solid surface, the segment den- 
sity profile and Green function only depend on the coordinate z so that 0(r) =(t){z) 
and G(r, ra) = G{z,n). In the (isotropic) bulk region G{z,n) Gb{n), i.e. inde- 
pendent of z. By solving the Edwards equation ()2.1|) with the initial condition 
in Eq. (j2.2j) . one may show that in the bulk region: Gt,{n) = ^-Ubn/kj^T ^ where we 
have defined Ub = U{z = 00). It is convenient to redefine the statistical weight to 
absorb this trivial n-dependence: 

Z{z,n) = e^''"/'=«^G(2,n). (2.5) 

The Edwards equation fEq. ()2.1|) ) and segment density fEq. ()2.3|) ) are then given 
by 



Z{z,n), (2.6) 



TV 



<f){z) = ^ JdnZ{z,n)Z{z,N-n), (2.7) 


with the initial condition: 

lim Z{z,n) = 1. (2.8) 

n— >0 

In the self-consistent field model, the external potential U{z) is in turn expressed 
in terms of the polymer segment density with the result that the set of equations 
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(|2.6|) - (j2.7|) becomes self-consistently closed. Various forms for the self-consistent 
external potential in terms of the segment density may be proposed. It should in 
some way describe the interaction between different polymer segments and of the 
polymer segments with the solid wall. A convenient form would be: 

The first two terms describe the polymer segment interaction in a virial expansion 
where vq, the so-called excluded volume parameter, is proportional to the second 
virial coefficient and Wq is proportional to the third virial coefficient of segment- 
segment interactions. The third term describes the non-local character of the 
interaction and reflects the non-homogeneous nature of the polymer segment dis- 
tribution. Such a term is analogous to the squared-gradient term in the van der 
Waals theory for the liquid-vapor interface. The last term in Eq. (|2.9p gives the 
interaction with the wall. In the following we will assume that it is short-ranged 
and that it can be approximated by a delta function located at the wall [1]: 

The parameter (i> is termed the extrapolation length [1]; its inverse is a measure 
of the surface interaction strength and leads to enhanced polymer adsorption, 
0(0) >(j)b- In the following, we will investigate the situation where 1/dis small [21] 
leading to weak adsorption. 

Weak adsorption 

The assumption of weak polymer adsorption implies that 6(f){z) = (f){z) — (pb'^'Pb 
and 5Z{z, n) = Z{z,n) — 1^1. Linearization of Eqs. ()2.6p and ()2.7|) then leads to: 

|-.Z(.„) - (2.n, 

b<\){z) = ^ [dn5Z{z,n), (2.12) 



N 
b 

with the initial condition: 

lim SZ{z,n) = 0. (2.13) 
The general form for the external potential expanded around (f) = (f)b is given by 
^ = vS^iz) -m5<P"iz) - ^ ^ Siz) . (2.14) 



The dimensionless parameter f is a generalized excluded volume parameter. Since 
we have expanded around = 06 rather than around = 0, as is done in the virial 
expansion in Eq. ()2.9|) . we have that v = v{(f)b) is not necessarily equal to the second 
virial coefficient {v ^ vq) and thus may contain higher order interactions. For 
stability, we do require f >0; to describe the situation t><0, it would be necessary 
to include more terms in the expansion around 0;,. The above expression for the 
external potential is used when we make the comparison with the SF calculations. 

For the situation that m = 0, it was shown [21] that the linearized set of 
equations in Eqs. (|2.1ip - (|2.14j) may be reformulated into a free energy functional 
formalism with the free energy functional given by 



a^F[5(t)] 
AkuT 



CO oo 



80. AT 



dz dz' S(j)'{z) 5(f)'{z') 



,z — z 
~2R7 



.z + z\ 
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+||d.[^0(.)^]-^l0(O). 



(2.15) 



where the polymer's radius of gyration -Rg = ^/WoFJq and where 



dk 
2^ 



ikx 



fc2 -4 + 4e-'=V4 • 



(2.16) 



Minimization of the free energy in Eq. ()2.15j) gives for the segment density profile 
[21] 



50(z) 



oo 

dt cos(-^) 



- 1 +exp(-t2 







(e/4) t^ + f^-l + exp(-t2) 



(2.17) 



where we have introduced the parameter e as the (square of) the bulk correlation 



length, ^b = a/\^3v (j)b, divided by the polymer's radius of gyration: 



The expression for the free energy can be generalized to the situation m^O 
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(2.19) 



The resulting segment density profile is then given by 



oo 



•'hit /"j. 



5<j){z) = dt cos( — ) 



TidRg J Rg 



t2 - 1 + exp(-t2) 







(e/4) + (1 + pt^) (t2 - 1 + exp(-t2)) 

(2.20) 

where we have introduced the parameter P as the (square of) the length defined 
by the ratio of m and v, ^m = \/ m/v, divided by the polymer's radius of gyration: 

The length is the typical length scale connected to spatial inhomogeneities 
in the segment density caused by the interactions between segments excluding 
interactions due to the chain's connectivity (which are described by ^b). In a good 
solvent C,b dominates over C,m while gains in significance near the theta region. 

The main goal in this article is to compare the segment density profiles given 
in Eqs. fl2.17|) and ()2.20j) to numerical results from the Scheutjens and Fleer lattice 
model. 



3 Scheutjens and Fleer self-consistent field model 

The Scheutj ens-Fleer model [18,19] gives an efficient way to solve, self-consistently, 
the discrete version of the Edwards equation ()2.(i|l on a lattice. The lattice dis- 
tance is taken to be equal to the polymer segment length a, and 1/A is the lattice 
coordination number that we set equal to 1/A = 6 corresponding to a cubic lattice. 
The distance z away from the surface is here a discrete variable, z = 0, 1, . . . , M 
layers (See Figure 1). A hard surface is located at z = 0, and at 2; = M we have 
reflective boundary conditions (M should be chosen large enough for the system 
to reach bulk density). 

The polymer is modeled as a one- dimensional walk perpendicular to the sur- 
face, where the walk can either move to the neighboring layer or stay in the layer 
it already is. Each step is weighted with the probability of going to a neighboring 
layer (A) or staying in the same layer (1 — 2A). This gives the following recurrence 
relation for the Green function G{z,s) (the discrete analog of Z{z,n)), which is 
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the statistical weight of finding the end segment s of a chain in layer 

G{z,s + 1) = G{z) [XG{z-l,s) + {l-2X)G{z,s) + XG{z + l,s)] , 

= G{z) <G{z,s)> . (3.1) 

As a starting point for the recurrence relation we have that G{z, 1) — G{z), the 
statistical weight of a single polymer segment. The brackets < . . . > denote the 
neighbor- weighted average as defined above. In the first layer, z = l, the neighbor- 
weighted average is defined as: 

<G'(l,s)>= (l-2A)G'(l,s) + AG'(2,s). (3.2) 

In analogy with the continuous case, the segment density profile is obtained 
by summing over all possible paths of a polymer of length passing through 
position z: 

<t>,^ G{z,s)G{z,N-s + l) 
<t>i^) -^1^ G{^^ > (3-3) 

.5 = 1 ^ ' 

where the factor G{z) in the denominator accounts for the double counting of 
the statistical weight at layer z. The statistical weight G{z) is written as the 
Boltzmann factor of the segment potential 

G{z) = c-^sfWAbT _ (3 4) 

The segment potential Usf{z) accounts for the interactions between polymer seg- 
ments and the surface-monomer interaction. The zero of the potential is at 2; = oo 
so that Usf{z — oo) — and G{z = oo) — 1. In the Scheutjens and Fleer model 

the form for the segment potential is chosen such that it is consistent with the 
Flory-Huggins free energy expression for a homogeneous polymer solution [13]. It 
is given by [18]: 

^ = - In (^T^) - X « m > - < 1 - m > +1 - 20.) - Xs 5.1 . 

(3.5) 

The first term is derived from the translational entropy of the "solvent molecules" 
(the sites not occupied by the polymer segments) in the Flory-Huggins theory. 

The second term describes the effective monomer-monomer interaction energy 
through the Flory parameter x [13]. (Notice that this term is non-local since the 



interactions between monomers in different layers is explicitly taken into account 
by the use of < (f){z) > instead of (f){z).) The final term describes the surface- 
monomer interaction through the Silberberg parameter Xs [26]. This parameter 
corresponds to the energy gain for a polymer segment to replace a solvent molecule 
at the surface {z = l). 

As input parameters for the SF calculations we thus have: 0^, A^, cind Xs- 
Many extensions to this model are possible and many have been made [2,20]. 
One may, for example, introduce different polymer architectures, inhomogeneities 
in more than one dimension, and curved lattice geometries, but the underlying 
principles remain the same. 

4 Comparison of results 

To make a comparison between the density profiles from the analytical theory and 
the Scheutjens-Fleer self-consistent field model, one needs to relate the parameters 
in the analytical theory to those in the lattice model. In particular, one needs to: 
(z) relate the interaction parameters v and m to the Flory-Huggins parameter x, 
and (a) relate the extrapolation length d (or rather 1/d) to the surface interaction 
parameter Xs- 

The discrete self-consistent field equations in the SF model can be transformed 
into continuous equations by assuming that the variations with z are small. In 
Appendix A we show that the recurrence relation in Eq. ()3.H) then reduces to the 
continuous Edwards equation ()2.6p with the external potential equal to Usf{z)- 
To make the correspondence between x and the parameters v and m, one therefore 
needs to expand the expression for Usf{z) in Eq. ()3.5p around (j) = (f)b, using that 
<0>— >0 + 0" in the continuous limit (cf. Appendix A), and compare the 

result with 6U{z) in Eq. ()2.14|) . One then finds that 



It is less straightforward to arrive at a relation between the extrapolation 
length d and the surface interaction parameter Xs- It should be realized that the 
presence of a solid surface is treated fundamentally different in the two models. 



l-0b 



1 



-2x, 




(4.1) 
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In the weak adsorption model, the treatment of the presence of the sohd surface 
is derived from the de Gennes model for polymer adsorption [1]. In the de Gennes 
model, the free energy functional is strictly defined for z>0; an interaction energy 
is added located at 2; = fEq. ()2.10|) ). but such a term merely results in a boundary 
condition to the differential equation obtained from minimizing the free energy 
functional. As a result, the polymer segment density does not necessarily go to 
zero at the wall, 0(0)7^0. 

In the SF model, it is implicitly assumed that the polymer segment density at 
the wall is zero. This can be read off from the definition of the neighbor-weighted 
average in the first layer (see Eq. ()3.2|) ). which for the segment density reads: 



i.e. 0(0) = 0. It turns out that a consequence of this condition is that enhanced 
adsorption only occurs in the SF model above some threshold value, Xs >Xsc- The 
situation Xs = Xsc then corresponds to the condition l/d = 0, which gives (f){z) = (j)f, 
for all z in the de Gennes model. The threshold value as a function of x 
therefore be determined by using that at the threshold one should have (j){z)^(f)b 
for all z. For very long chains one finds: Xsc = — Ax — ln(l — A) [27]. Next, 
by comparing the expressions for the wall interaction terms in the two models, 
Eqs. (j2.14|) and (j3.5|) . one arrives at the following identification [9,27] 



where we used that i< — >a (1 — A) 6{z) [27]. 

We now have the necessary ingredients to compare the analytic expressions 
for the segment density profiles in Eqs. ()2.17|) and ()2.20|) to the numerical results 
of the SF model. It is good to realize that the inverse of the extrapolation length, 
1/d, only shows up as a prefactor to the analytic expressions for the segment 
density profile. Therefore, if one is interested in only the qualitative features of 
the segment density profile, knowledge of the precise value of l/ci is not needed. 
Only when one wants to make a quantitative comparison with the SF model, it is 
necessary to use the above relation between 1/d and Xs- 

Next, we first consider the case of an athermal solvent (x = 0) and compare 
the SF model results to Eq. (|2.17|) . In particular, we will investigate the distal 



<0(1)>= (1-2A)0(1) + A0(2), 



(4.2) 



1 





(4.3) 
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oscillations for weak and strong polymer adsorption. Second, we consider a tlieta 
solvent (x= 1/2), as was previously done by van der Guclit et al. [22], who showed 
that when the bulk density is below some threshold value, the oscillations in the 
distal profile disappear similar to the Fisher- Widom transition [24, 25] in the 
damped oscillations for a simple liquid near a surface. A quantitative comparison 
is made with the analytical expression for the segment density profile in Eq. ()2.20|) . 
Finally, we consider the location of the Fisher-Widom transition with respect to 
the complete Flory-Huggins phase diagram [13] and also make a comparison with 
MC simulations [15]. 

4.1 At hernial solvent 

For an athermal solvent we have that x = so that: 

^ 1 , m = . (4.4) 



In Figure 2a, the polymer segment density profile as calculated by the SF model 
is shown for two polymer chain lengths (A^= 1000 and A^ = 10000) in the case of 
weak polymer adsorption (xs = 0.20). To show 6(f){z)/(f)b on a logarithmic scale, 
here, and in later plots, the absolute value has been taken. Also shown as the 
solid line is the ground state dominance segment density profile [1]: 

(4.5) 



with the integration constant Xq = (1/2) arcsinh(2(i/^fe). In the comparison of 
Figure 2a, the expression for a/d in Eq. (j4.3p was used. 

The profiles corresponding to the two chain lengths are close to each other 
and close to the ground state dominance profile in this distance range. At very 
large distances, however, markedly different behavior can be observed; in Figure 
2b, the same SF model results as in Figure 2a are plotted, showing oscillations 
in the segment density profile on the order of the polymer's radius of gyration. 
The two solid lines in Figure 2b are the density profiles calculated by Eq. (|2.17|) . 
using Eq. (|4.3|) . Without any adjustable parameters, a strikingly good quantitative 
agreement between the SF model and Eq. ()2.17|) is obtained. 
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Even for strong adsorption the density profile in Eq. (j2.17p is expected to give 
an accurate description for distances far away from the surface, since also there we 
have that 6(f){z) <^(j)b. However, since the density profile differs at short distances, 
the boundary condition may need to be replaced by some effective boundary 
condition. The result is that although the shape of the density profile is correctly 
described by Eq. fl2.17|) for strong adsorption, it may need to be shifted in the 
^-direction introducing the shift as an additional fit parameter. In Figure 3, the 
polymer segment density profile as calculated by the SF model is shown for a long 
polymer chain (A^ = 1000) and for two values of the surface interaction parameter, 
Xs = 0.25 and Xs = 1? corresponding to strong adsorption. The solid lines are 
the density profiles calculated by Eq. ()2.17|) . using Eq. ()4.3|) to set the scale of the 
y-axis. Furthermore, the density profiles are shifted in the 2;-direction by —1.5a 
(lower solid line) and —4.5a (top solid line). The asymptotic profiles are seen to 
indeed match quite well. 

4.2 Theta solvent 

We next consider the case of a B-solvent (x = l/2) where we have that 

1 

v = - — -1^0b, ^=-^- (4-6) 

1-06 6 

It is good to realize that theta solvent conditions do not correspond to a vanishing 
generalized excluded volume parameter (v^O). 

In Ref. [22] the observation was made that the oscillations, prominently present 
in Figures 2 and 3, disappear in a G-solvent when the density is below a certain 
threshold bulk density, similar to the Fisher-Widom transition [24] for a simple 
liquid. In Figure 4, we show the characteristic signature of the Fisher-Widom 
(FW) transition going from 0b = 0.01 (solid line) to (pb = 0.001 (dashed line), 
for = 1000 and for the same surface interaction strength, Xs = —1/12, as in 
Ref. [22]. (Note that this value is below the adsorption threshold value, so that 
we are dealing with polymer depletion instead of adsorption.) To understand 
this qualitative change in the behavior of the segment density profile, we further 
investigate the expression for the density profile in Eq. (j2.20j) . 
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In Appendix B, it is shown that the asymptotic behavior is determined by the 
poles in the complex plane of the integrand's denominator of the expression for the 
density profile in Eq. ()2.20|) . The resulting expression for the asymptotic segment 
density profile as given by Eq. ()B.6|l shows that the relevant poles correspond to 
either a purely exponentially decaying function: 

S^iz) oc exp(-^) , (4.7) 

or to an exponentially decaying sinusoid: 

S<P{z) (X exp(-^) sin(§^ + ^) . (4.8) 
Kg Kg 

Which of these two functions dominates the asymptotic behavior depends on the 
relative magnitudes of Aq and Ai; when Aq > Ai, the segment density decays 
asymptotically as an exponentially decaying sinusoid whereas it decays exponen- 
tially when Aq<Ai. The value of the coefficients Aq, Ai and Bi, depend on the 
value of the parameters e and (3. In Figure 5a, we have set (3= 0.18 and plotted 
Aq, Ai and Bi as a function of e. At a certain value of e (in this example: ^ 
6.464), Aq and Ai cross marking the location of the Fisher- Widom transition 
(arrow). By varying (3, the whole locus of Fisher- Widom transitions can then be 
traced in this way (see Figure 5b). 

To compare the analytically determined Fisher- Widom transition in Figure 5b 
with the results from the SF model, we located the Fisher-Widom transition by 
determining the bulk density where the first oscillation in the numerical density 
profiles disappears. This is the point where the first minimum (located closest to 
the wall) of the oscillating density profile crosses the value of the bulk density. 
This procedure gives only an approximation to the location of the real Fisher- 
Widom transition. A closer investigation reveals that the density oscillations do 
not disappear all at once at the same polymer bulk density; the first oscillation 
disappears before the second oscillation with decreasing bulk density, and presum- 
ably before all oscillations beyond that. Still, since the bulk densities at which 
subsequent oscillations vanish differ only slightly, we believe this procedure to give 
a fairly accurate approximation. 

In Figure 6, the polymer segment density profile is shown for a 0-solvent with 
= 10000 for various bulk densities just below and above the Fisher-Widom 
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transition: 0^ = 0.00078 (solid line), (pb = 0.00082 (dashed line), (pb = 0.00086 
(dot-dashed line), and (f)b = 0.00090 (dotted line). In (a), the density profiles are 
calculated in the SF model; in (b) the density profile is calculated using Eq. ()B.6|) . 
Good quantitative agreement is obtained. 

In Figure 7, we have located the Fisher-Widom transition in the SF model 
by determining the bulk density where the first oscillation disappears for different 
chain lengths (circles). The surface interaction strength Xs = as in Ref. [22]. 

The numerical data are consistent with a scaling of the bulk density (pl ~ 
The solid line is the Fisher-Widom transition shown in Figure 5b, where we have 
used that = e/{2[3'^) and (pb = 2/3/e to transform from (e, (3) to (A^, (pb) as 
variables. Good agreement is obtained with the SF calculations. 

Also shown in Figure 7, as the squares, are the SF calculations results for A^ = 
200 and A^ = 1000 reported in Ref. [22], which differ somewhat from our results. 
Since these SF calculations are done with exactly the same set of parameters used 
to determine our results represented by circles, this difference must be attributed 
to the different procedure used in Ref. [22] to locate the Fisher-Widom transition. 
In the analysis in Ref. [22] the segment density profile is taken to be of the following 
form: 

— - — = tanh (- ) - 1 + Cnum exp( — ) sm( — h v^num) • (4.9) 

Then, by locating, as a function of 06, the position where Cb,num/2 equals i?g/Anum, 
the Fisher-Widom transition is determined. 

Even though the analysis in Ref. [22] is done in a more judicious manner than 
simply fitting the numerical segment density profile to Eq. ()4.9j) (see Ref. [22] for 
details), the large number of parameters involved makes it difficult to accurately 
determine the location of the Fisher-Widom transition — especially since at the 
transition the relevant length scales cross. For instance, in the example shown in 
Figure 4, the Fisher-Widom transition as determined by the bulk density where 
the first oscillation disappears, is located at 0^^ ~ 0.00707, whereas the value 
reported in Ref. [22] is 0^^ ~ 0.032. This result would imply that the density 
profile depicted by the solid line in Figure 4 for (pb = 0.01 lies in the purely 
exponentially decaying-region, which seems unlikely. 
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Nevertheless, it is expected that this fitting procedure is more rehable in the 
region where the different length scales are well-separated, e.g. when 0^^^^^. 
In this region, the fit parameters were numerically determined as [22] 

1 P 27r 0.52 a 

'^""^ 01971 '^"'"^ 1:571^ ^ ' 

These values can be compared to the analytical density profile in the corresponding 

regime: e and (3 <^1. For (3 = and expanding in £, it was shown that the 

segment density profile is then given by [21] 

.|e--/^. + Ce.p(-^) =in(|^ + »=), (4.11) 

with 

A = 2.217792... , 5 = 1.682188..., = — ^ ~ ^^4^ • (4-12) 



>b 



Thus, we find in the region (ph ^ good agreement between the analytical 
results [21] and the fitting parameters in Eq. ()4.10|) . 



4.3 Phase diagram 

Since the Fisher-Widom transition found in a 0-solvent seems to be absent in 
an athermal solvent, it is worthwhile to locate the Fisher-Widom transition as a 
function of solvent quality. 

In Figure 8, the Fisher-Widom transition is shown with respect to the full 
polymer phase diagram for = 100. As a reference, the Flory-Huggins spinodal 
and binodal regions (dashed lines) are also shown [13]. The symbols are the SF 
results for the Fisher-Widom transition for enhanced polymer adsorption Xs = 0.5 
(circles) and for polymer depletion Xs = (squares). The solid line is the Fisher- 
Widom transition plotted in Figure 5b transformed from {e, (3) to (x, (t>b) as 
variables, using Eq. ()4.H) . Good agreement is obtained with the SF calculations, 
independent of the polymer-wall interaction strength. 

In the derivation of the analytical result for the Fisher-Widom transition, we 
have assumed that the generalized excluded volume parameter is always positive, 
f >0. To verify that the Fisher-Widom transition lies within the region of appli- 
cability, the locus of points where f = (which gives 2%= 1/(1 — 0^)) is shown in 
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Figure 8 as the dotted line. Furthermore, the approximate location of the dilute 
to semi-dilute transition {(f)* = 1/ (vN) [28], with v given by Eg . (14.11) ) is shown as 
the dot-dashed line in Figure 8. Since both the mean-field model presented here 
as well as the SF model are expected to be reliable only in semi-dilute and concen- 
trated solutions, the location of the Fisher-Widom transition shown is physically 
most relevant in the theta region. 

As a final remark, we mention that a useful approximation to the location of 
the Fisher-Widom transition can be obtained by approximating by its e 
value: ^ o.20 (see Figure 5b). This gives 2x^'^ ~ N/[{N + 5) (1 - (p^)]- For 
very long polymer chains, the Fisher-Widom transition thus approaches the line 
where f = 0. 

4.4 Comparison with Monte Carlo simulations 

Recently, de Joannis et al. [15] performed Monte Carlo (MC) simulations of poly- 
mer chains at surfaces in an athermal solvent. In Figure 9a, the reduced segment 
density 6(f){z)/(f)b is shown as a function oi z/Rg (Figure 9b shows the same profiles 
on a logarithmic scale). For the MC simulations (circles), the following parameter 
values were used: chain length A^ = 200, bulk density 0?, = 0.0216, radius of gyra- 
tion i?g/a = 9.76, and surface interaction strength es= 1.0 k-oT [15]. The squares 
are the results from the SF model for the same parameters (A^ = 200, x = 0) bulk 
density 0^ = 0.0216) taking Xs = l-0 (variation of Xs does not lead to very different 
results). Shown as the solid and dashed line is the density profile calculated using 
Eq. ()2.17p for a/c/= 7.35 and a/(i = 4.09, respectively (both profiles are shifted in 
the 2;-direction to fit the data). The value of a/d= 7.35 was chosen such that 
the depth of the minimum in the segment density profile matches that of the MC 
simulations, whereas the value of a/(i = 4.09 was chosen to correspond to Xs = l-0 
(using Eq. ()4.3|) ). The SF model corresponds well with the MC results although 
the depth of the minimum differs by a factor of two. The asymptotic formula in 
Eq. ()2.17p is able to describe the MC-data (two fitting parameters) and SF model 
(one fitting parameter) well. So, also in the Monte Carlo simulations oscillations 
in the polymer concentration profile can be observed, even though it is not pos- 
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sible from the numerical data (and it is hard to get more accurate profiles) to 
observe more than the first oscillation. 

5 Discussion 

In this article, we have compared the analytic expressions for the segment density 
profiles obtained in the context of the weak adsorption model [21] to the numer- 
ical results of the Scheutjens-Fleer model [18,19]. The analytic expressions are 
expected to be valid when the polymer adsorbs only weakly to the solid surface 
or, for strong adsorption, when the distance to the surface is large. Indeed, excel- 
lent quantitative agreement for the distal, oscillatory density profile was obtained 
for very good (athermal) solvent conditions, for both weak and strong polymer 
adsorption (Figures 2 and 3). 

We have further investigated the previously observed Fisher-Widom transition 
in SF model calculations by van der Gucht et al. [22] for a G-solvent. We demon- 
strated that the location of the Fisher-Widom transition can be described by the 
weak adsorption model when the non-local character of the monomer-solvent in- 
teraction is taken into account. The quantitative agreement that is then achieved 
(see Figures 6, 7 and 8), gives confidence in identifying the non-local character 
of the interactions in the SF model as the source of the Fisher-Widom transi- 
tion observed [22]. To further substantiate this identification, we also carried out 
separate SF calculations in which the non-local interactions were removed by re- 
placing <(j){z) > by (j){z) in Eq. ()3.5|) . In this way, we could explicitly verify [29] 
that without the non-local interactions present, the oscillations always remain 
with decreasing bulk polymer concentration while maintaining the quantitative 
agreement with the segment density profiles of the weak adsorption model. 

Although it is shown that the introduction of a non-local interaction term is 
useful for the comparison with the SF model, one should address the relevance 
of it when comparing with real, experimental polymer systems. Both the model 
presented here as well as the SF model are based on the mean-field assumption 
which is expected to be reasonable for semi-dilute and concentrated solutions 
(overlapping polymer coils) but is expected to fail for dilute solutions (isolated 
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coils), where fluctuations and inhomogeneities are more pronounced. Under good 
solvent conditions the Fisher-Widom transition occurs at polymer concentrations 
well below the dilute to semi-dilute transition, ^ 0* ~ 1/{vqN) [28], and 
the mean-field results for the location of the Fisher-Widom transition are prob- 
ably not so reliable. In theta and poor solvent conditions, close to the spinodal 
decomposition region (theta region) , the Fisher-Widom transition shifts to much 
larger concentrations and occurs at a concentration comparable to the dilute to 
semi-dilute transition concentration, ~ 0* ~ 1/(^0-^"'^^^) [28] (see Figure 8). 
We therefore expect our mean-field results to be most relevant to describe the 
occurrence and location of the Fisher-Widom transition in experimental polymer 
systems in the theta region. 

It may seem somewhat surprising that the non-local interactions in the SF 
model, whose range is only one lattice distance a, may interfere with the oscillatory 
behavior on the scale of the polymer's radius of gyration. However, the length 
scale connected with the non-local character is rather = \frnfv which diverges 
(similar to the bulk correlation length near the critical point in an ordinary liquid) 
when V is small, i.e. when the solvent is poor enough. 

The Fishcr-Widom transition in a simple liquid occurs when the attractive 
part of the interaction potential between the molecules dominates over the re- 
pulsive part [24, 25] . For high densities the hard core interaction gives damped 
oscillations away from the surface, while these oscillations disappear when the 
inter-molecular attraction balances the repulsion. The similarity with the transi- 
tion observed for polymers is striking; the 0-point defines the system conditions 
where the monomer attraction cancels the monomer repulsion. For a good solvent 
the effective monomer interaction is repulsive and oscillations are predicted for all 
densities [21]. There are, however, nontrivial differences between monomers and 
liquid molecules; the most important difference, of course, is the polymer chain 
connectivity. The chain connectivity gives rise to oscillations with a wavelength 
of the size of polymer coil, while in a simple liquid the wavelength is given by the 
molecular radius. It is tempting to think of the polymer solution as a "liquid" 
of polymer coils close to the surface. But this intuitive idea fails to consider the 
strong interdigitation of the coils at semi-dilute concentrations and that the coil 
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changes its conformation and size close to the surface (while a liquid molecule 
keeps its hard core shape) . One should therefore be careful to make too much out 
of such analogies. 

In the comparison of the Scheutjens and Fleer model [18, 19] to analytical 
theories such as the de Gennes model [1,4], and extensions thereof such as the 
weak adsorption model [21], a number of differences can be distinguished. A 
clear difference is the fact that the SF model is a lattice model with the obvious 
consequence that the parameter describing the distance to the substrate is discrete; 
^ = 1,2,3,.... The discreteness of space is most consequential in the very vicinity 
of the wall, where the monomer length scale becomes important. 

Another important difference between the SF model and the de Gennes-hke 
models concerns the boundary condition at the wall. In these latter models, the 
free energy functional is only strictly defined for z > 0. An interaction energy 
is defined located at ^ = 0, but the addition of such a term merely serves as a 
boundary condition to the differential equation which results from the minimiza- 
tion of the free energy functional. The result is that the polymer segment density 
not necessarily goes to zero at the wall, 0(0) 7^ 0. Even though in real poly- 
mer systems, the polymer segment density should become zero (or exponentially 
close to it) at some point, it is expected that a microscopic shift in the precise 
location of the wall — this location is not well-defined within a microscopic dis- 
tance anyway — should restore the vahdity of the de Gennes-like models at least 
for enhanced polymer adsorption [4]. 

In the SF model, the free energy functional is already defined for whole space 
from the start. The reason that the density profile is zero when 2; < is due to the 
explicit addition of an external potential which gives an infinite energy penalty for 
2; < 0. Apart from this external potential there may be an energy contribution to 
the first layer at z = l which will generally be attractive. The result of the infinite 
repulsion at ;s < is that the polymer density is zero at the wall 0(0) = 0. The 
presence of this extra boundary condition makes the SF model essentially different 
from the de Gennes-like models, but similar to the single chain mean-field theory 
by Szleifer and others [30,31] and computer simulations [14-16]. This boundary 
condition (but without the presence of any additional attraction) is also used in 
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the calculations by Eisenriegler and coworkers [3,32,33]. Still, it seems, as we 
have done here, that a meaningful comparison between the SF model and the de 
Gennes-like models can still be made when the adsorption energy in the first layer 
is above some given threshold value. 

A Discrete and continuous Edwards equation 

In this appendix, we show [2] how the recurrence relation for G{z, s) in Eq. p.ip 
in continuous form reduces to the Edwards equation ()2.6|) . First, we rewrite the 
recurrence relation for G{z, s) as: 

G{z, s + 1) e^sF(-)ABT _ Q^^^ s^^x[G{z-l,s)-2 G{z, s) + G{z + 1, s) ] , 

(A.l) 

where we used the expression for G{z) in Eq. ()3.4|) . To go from the discrete descrip- 
tion to the continuous one, we assume that G{z, s) is a slowly varying function of 
z and s. This means that: 

G{z, s + l) = G{z, s) + -^ G{z, s) + . . . , (A.2) 

and, keeping in mind that the distance z was rescaled by a, that 

G{z±l,s) G{z±a,s) =G{z,s)±a-^G{z,s) + ^-^G{z,s) + ... . (A.3) 

When G{z, s) is a slowly varying function, it implies that the potential Usf{z) is 
close to its asymptotic value Usf{z oo) = 0, so that we can expand 

Combining Eqs. ()A.2j) - ()A.4|l with Eq. ()A.l|l then gives: 

^ s) = Xa'^ s) - s) , (A.5) 

which, for A = 1/6, is the Edwards equation in Eq. ()2.6|) . 
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B Asymptotic behavior of the density profile 



In this appendix, we investigate the asymptotic behavior of the polymer segment 
density profile given in Eg. ()2.20| 1. We write Eg. ()2.20|) as: 



5<P{z) 



' dt cos{xt) Q{t) 



71 dR. 



(B.l) 



(B.2) 



(B.3) 



where we have defined x = z/Rq and Q(t) = f(t)/g(t) with 

fit) ^ t'-l + e-'\ 

git) ^ ^ + (l + /3t2) (t^-i + e-*^). 

The integral in Eg. ()B.l|) can be solved [34] 

oo 

Jdt cos(xt) Q{t) = vri ^ R' , 



where 

^ R' = sum of residues of Qiz) e*^^ in the upper half plane. 



The poles are determined by the zero's of g{z) in the complex plane (see also [25]). 
It turns out that there are an infinite number of poles in the upper half plane. 
Luckily, to describe the asymptotic behavior of Eg. ()B.l|) . it is sufficient to locate 
the pole with the lowest imaginary part. Depending on the value of e and (3, this 
pole is either purely imaginary — let's denote it as Zo = iAQ — or it is complex and 
then it comes in pairs: Zi = ±Bi + iAi, where Ao, Ai, and Bi are positive, real 
numbers. Taking only these three competing poles into account, we can determine 
the residues with the result: 



^ R' = z Co e-^» ^ + 2i Ci e-^i ^ sin(Si x + ^) + 



where 



Cn 



fi^o) 



fiZl) 



9'i^i) 



(p = arg 



fi^i) 



g'izo) ' 

and where zi is taken to be the pole in the upper right guadrant. 
As a final result, we thus have for the asymptotic density profile: 



(B.4) 



(B.5) 



Co , Aqz^ Axz . Biz 

" " + Cl exp( — sm( 



exp( 



Rp 



Rp 



^) 



(B.6) 
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Figure 1: An adsorbed polymer on a lattice. An impenetrable wall is located at 
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(b) 

Figure 2: Polymer segment density profiles 6(f){z)/(f)b for an athermal solvent (x = 
0) and for (j)b = 0.01. The symbols are the SF calculations for A^= 1000 (circles) 
and A^ = 10000 (squares) with Xs = 0.20. In (a) the solid line corresponds to the de 
Gennes profile Eq. fl4.5p . In (b) the solid lines are the density profiles calculated 
with Ea. (|^T7jl . 
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Figure 3: Polymer segment density profiles 6(f){z)/(f)h for an athermal solvent (x = 
0) and for (f)b = 0.01. The symbols are the SF calculations for Xs = 0.25 (circles) 
and Xs = 1-0 (squares) with = 1000. The solid lines are the density profiles 
calculated with Eq. (j2.17|) and shifted horizontally by —1.5 a (bottom solid line, 
corresponding to Xs = 0.25) and —4.5 a (top sohd line, corresponding to Xs= 1-0). 
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Figure 4: Polymer segment density profiles S(j){z)/(j)b for a 0-solvent {x = 0.5) 
calculated in the SF model; Xs = — 1/12 and A^= 1000. The bulk density (f)b = 0.01 
(solid line) and (j)b= 0.001 (dashed line). 



27 



10'^ lO' 10° lO' 10" 10"' lO* 

e 

(a) 




-3 -'' -1 I 3 4 

10 10' 10 10 10 10" 10 10 

e 

(b) 

Figure 5: (a) shows the coefficients Aq (dashed curve), and Ai, Bi (sohd curves), 
that determine the exponential decay and the oscillation period of the polymer 
density profile (see Eqs. ()4.7p and ()4.8|) ). as a function of e = 2/{v(f)bN) for /? = 
Qm/{Na'^v) = 0.18. At 5fw = 6.46398 . . ., and Ai cross. In (b) the locus of 
Fisher- Widom transitions is shown as a function of e and f3. 
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Figure 6: Polymer segment density profile 5(l){z) / (pi, for a 6-solvent with A^= 10000 
for various bulk densities: 0;,= 0.00078 (solid line), 0;,= 0.00082 (dashed line), 
0^= 0.00086 (dot-dashed line), and 0b = 0.00090 (dotted line). In (a), the density 
profiles are calculated in the SF model, using Xs = —1/12; in (b) the density 
profiles are calculated using Eq. ()B.6|) . 
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Figure 7: Bulk density corresponding to the Fisher- Widom transition for a ©- 
solvent {x = 0-5) as a function of polymer chain length A^. Symbols are the SF 
calculations for Xs = — 1/12; circles correspond to the present calculations, squares 
are the results reported in Ref. [22]. The sohd hue is the analytical result (see 
text for details). 




Figure 8: Fisher-Widom transition in the polymer phase diagram for N = 100. 
The symbols are the SF calculations for = 0.5 (circles) and Xs = (squares). 
The solid line is the analytical result (see text for details). The dotted line is 
where the generalized excluded volume parameter v — 0. The dilute to semi-dilute 
transition is shown as the dot-dashed line {(/)* ^ l/{vN)). Also shown are the 
Flory-Huggins spinodal and binodal regions (dashed lines). 
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(b) 

Figure 9: Reduced polymer segment density profiles (A^ = 200, (f)b = 0.0216). 
Circles are MC results [15]; squares are the SF calculations for Xs = 1-0. Solid line 
is the density profile calculated by Eq. ()2.17|) for a/d = 7.35 and shifted by —2a to 
fit the MC results. Dashed line is the density profile calculated by Eq. ()2.17p for 
a/(i = 4.09 (corresponding to Xs = l-0) and shifted by —3.5 a to fit the SF results, 
(b) shows the same results on a logarithmic scale. 
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